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ABSTRACT 

Context. V838 Mon is an eruptive variable, which exploded in 2002. It displayed the most spectacular light echo ever observed. 
However, neither the origin of the reflecting matter nor the nature of the 2002 outburst have been firmly constrained. 
Aims. We investigate the nature of the CO radio emission detected in the field of the light echo. In particular, we explore its connection 
to the echoing dust around V838 Mon. 

Methods. We observed the echo region in multiple CO rotational transitions. We present and analyse maps of the region obtained in 
the '^CO(l-O) and (3-2) lines. In addition, deep spectra at several positions were acquired in '-CO(l-O), (2-1), (3-2), and '^C0(1- 
0), (2-1). Radiative transfer modelling of line intensities is performed for chosen positions to constrain the kinetic temperatures and 
densities. We derive global parameters (e.g. mass, distance, total column density) of the emitting cloud. 

Results. We found that a compact molecular cloud is located within the echo region. The molecular emission is physically connected 
to the dusty environment seen in the optical echo and they both belong to the same translucent cloud. The interstellar nature of the 
cloud is confirmed by its high mass of 90-150 Mq. We propose that the cloud consists of material remaining after the formation of 
the cluster to which V838 Mon belongs. This indicates that the eruptive star has young age (3-10 Myr). 

Key words, radio lines: ISM - ISM: clouds - ISM: molecules - stars: individual: V838 Monocerotis 



1. Introduction 

V838 Mon is an unusual variable, wh ose nova-like eruptio n was 
observed at the beginning of 2002 (iMunari et al ] 120021) . The 
object displayed a peculiar light curve with at least three dis- 
tinct maxima seen in the optical. At maximum light, V838 Mon 
reached a luminosity of 10^ L0 and displayed the spectrum of 
an F-type star. The outburst lasted for about 3 months and was 
ended with a steep decline, e.g., a decline in V-band brightness 
of about 8 mag within one month. The decline was accompanied 
by a rapid decrease in the effective temperature (Tvlenda 2005) 
and the object eventually tu rned into a cool s upergiant with a 
spectral type later than MIO (lEvans et al.ll2003h . V838 Mon has 
retained its low tem perature since then and its spectral type i s 
now close to M 6-7 (iRamifiski e t alJl2009l : lTvlenda et aljr201 Ih . 
The observational characteristi cs of the object during and after 
the outburst can be found in lCorradi & Munaril(l2007b . 

Several different explanations of the eruption of V838 
Mo n have been proposed, which were critically discussed 
by iTvlenda & SokeJ 32006). The most promising model ap- 

J iears to be that of a stellar collision fo llowed by a merger 
Soker & Tvlendal2003HTvlenda & Soke32 006). A crucial point 
in establishing the nature of the outburst is to constrain the evo- 
lutionary status of the progenitor of V838 Mon. 

V838 Mon displayed a spectacular light echo, which has 
been monitored since 2002 February. Multi-epoch images of 
the echo show a diffuse neb ula surrounding the star that ex - 
hibits a complex morphology (iBond et al.l2003l:ICarlqvistl2005h . 
Polarimetric observations of the echo were u sed to derive a ge- 
ometric distance of 6 kpc (ISparks et al.ll2008l) . The natu re of the 
dusty material responsible for the echo remains unclear. iTvlendal 



(I2OO4 see also iTvlend a et al.ll2005h interpret it as an interstel- 
lar cloud, while iBondl (12007) advocates a circumstellar origin. If 
the echo material was ejected by the progenitor of V838 Mon, 
its properties (e.g. its mass ) should be consistent with the nature 
of the progenitor object. Constraining the nature of the echoing 
cloud is therefore important to more clearly understand the enig- 
matic eruption from 2002. 



Emission in the lowest CO rotational transitions was de- 
tected at the position of V838 Mon and in its close vicinity 
(Kamitiski et al. 2007; Deguchi et al. 2007). First observations 
of the molecular line emission suggested that the source is ex- 
tended. These observations were however of poor spatial resolu- 
tion or incomplete sky coverage, so the connection of the molec- 
ular gas to either V838 Mon or the light-echo material was un- 
clear. .Kamitiski (2008) observed 13 positions within the echo 
region with the IRAM 30-m telescope and found that the emis- 
sion is indeed extended and shows a complex spatial distribution. 
However, its relation to the light echo material remained unclear. 



In this paper, we present comprehensive, follow-up observa- 
tions of the echo region in several CO rotational transitions. The 
observations are described in Sect.|2]and the results are presented 
in Sect. |3] The observations reveal the presence of a molecular 
cloud within the echo region. We constrain its basic physical pa- 
rameters in Sect. |4] In Sect. |5] we show that the molecular cloud 
is physically related to the dusty medium seen in optical echo 
images. In Sect.|6] we discuss the nature of the cloud. 
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2. Observations and data reduction 

A restiicted region around V838 Mon was mapped in the 
'^CO(l-O) and (3-2) transitions, for which the observations, 
data reduction, and map assembly procedures are described in 
Secs. l2.ll and l2.3l In addition, some positions within the mapped 
region were observed in the standard on-off (position switch- 
ing, PSW) mode. These observations include integrations in the 
12CO(1-0), i^CO(l-O), and HCO+(1-0) Hnes. Technical details 
about these measurements are given in Sect. 12.21 

In this paper, w e often make a ref erence to observations pre- 
viously reported in iKamihs kil (l2008h . We summarize these data 
in Sect. 12.41 Some figures in this paper also include an archival 
optical image of the light echo performed with the HST The 
source of the image and the applied reduction procedure are 
shortly described in Sect l2.6l 

All velocities in this paper are given with respect to the local 
standard of rest (LSR) and celestial coordinates are given for the 
epoch J2000.0. 

2.1. BEARS map in ^^CO(1-0) 

Mapping in '^CO 7=1-0 ( 115.27GH z) was per formed with the 
BEARS multi-beam array (ISunada et al. 2000; Yamaguchi et al.l 
[2000) installed on the 45-m telescope of the Nobeyama Radio 
Observatory (NRO), Japan. The array consists of 25 detectors 
arranged in a 5 x 5 grid separated by 41'.' 1 from each other 
The half-power beam width (HPBW) of the 45-m telescope at 
1 15 GHz i s 14'.'9. The map wa s obtained in the on-the-fly (OTP) 
technique dSawada et al.ll200 8') with 5" spacing between succes- 
sive scan rows. The whole region of interest was observed sev- 
eral times and successive coverages were obtained in orthogo- 
nal directions. The position angle of the array during observa- 
tions was PA = 58° east of north. We mapped the vicinity of 
V838 Mon in two steps. First, in 2008 April 2^ we observed an 
area of 320" x 320" centerd at the position a = 07''04'"05';70, 
6 = -03°50'25'.'0. In this area, only the inner part of 137"x 137" 
was fully sampled. Next, on 2008 April 12 we mapped a larger 
area of 640" x 640" with full sampling within the inner region 
of 457" X 457". 

The central off position used for the OTF observations was at 
a = 07''05"'25!70, 6 = -04°03'25'.'0. This position was found to 
be free of CO(l-O) emission in the sensitive observations with 
the Delingha telescope repotted in Kamifiski et al. (2007). 

All the BEARS data were ac quired with the ar ray of digi- 
tal autocorrelation spectrometers (ISorai et aDl2000b in a mode 
that provides a bandwidth of 32MHz (83km s"') per detector 
and a spectral resolution of 31.25kHz (0.081 km s '). The sys- 
tem noise temperature, including atmospheric contributions, was 
typically T,y, = 400K (DSB). 

Data were calibrated wit h the standard chopper-wheel 
method (lKutner&Uhchlll981h . The calibration of the antenna 
temperature, T^, was controlled by observations of a strong CO 
source in the Orion KL region, located at a = 05''35'"14!5, 
6 = -05°22'30'.'4. Day-by-day measurements demonstrated that 
the calibration was stable to within 14% (3cr noise). Pointing, 
crucial for high-frequency observations at the 45-m telescope, 
was performed on strong SiO maser sources every ~1 h. The ob- 
servations were carried out in slow wind (i.e. with an average 
wind speed lower than 5 m s ') or no wind conditions and the 
pointing should be accurate to within ~ 9" (3cr). 

Data reduction and map g eneration were performed with 
NOSTAR dSawada et al.ll2008r . Baselines of up to third order 
were subtracted from the spectra and the data were rescaled to 



Table 1. NRO 45-m telescope specification. 



transition 


frequency 


HPBW 


'7mb 




[GHz] 


["] 




'-CO y=i-o 


115.271204 


14.9 


0.32 


"CO y=i-o 


110.201353 


14.9 


0.40 


HCO+ 7=1-0 


89.188518 


18.2 


0.44 



single side-band (SSB) scale using conversion factors provided 
by NRO. The data were then converted to main beam bright- 
ness temperature units, Tmb - T^/rjmb, using the main beam effi- 
ciency of 77,„ft=0.32. To reduce the scanning noise effect, the data 
were regiidde d and conv olved into a map using the basket-weave 
method ( Mangum et al.li2007 ). The data were convolved with a 
Gaussian-tapered first-order Bessel function (BesselxGauss) as 
the gridding convolution function. For a grid spacing of 7'.'5 in 
the final map, the described procedures resulted in an effective 
beam-size of 19'.'2 (HPBW). 

2.2. Single position measurements witli BEARS 

Several positions within the echo region were observed in '^CO 
7=1-0 (110.20GHz), HCO+ 7=1-0 (89.19GHz), and '^CO 
7=1-0 in the standard PSW mode. The measurements were per- 
formed with BEARS. The HPBWs and main-beam efficiencies 
for the NRO 45-m telescope at the observed frequencies are 
given in Table [T] 

Measurements in '^CO and '^CO were obtained for two cen- 
tral positions (note that the actual number of points observed 
with BEARS is 25 for every specified central position). One 
of the observed positions was the emission maximum found 
on the '^CO(3-2) map (see Sect. I2.3l l. which we hereafter call 
P32. It is located at a = 07''04'«05^.40, 6 = -03°50'13'.'0. 
We also observed the position with '-CO(l-O) maximum, as 
identified in the BEARS map. This position has coordinates 
a = 07''04'"05^;52, 6 = -03°50'00'.'0, which we refer to here- 
after as PIO. In the case of the HCO+(1-0) fine, only PIO was 
observed. Additional observing details for the PSW spectra are 
listed in Table |2] 

The calibration and data reduction of the PSW observations 
was completed using analogous procedures as for the OTF ob- 
servations with BEARS (see Sect. 12. 111. Data reduction w as per- 
formed using the NEWSTAR package dlkeda etalj|2001h . 

2.3. Map in ^^CO(3-2) 

The map in '^CO 7=3-2 (345.7956GHz) was obtained with 
HARP (Heterodyne AiTay Receiver P rogram for the B-band) 
dBuckle etalJ 120091: ISmithet al.ll2008h newly installed on the 
James Clerk Maxwell Telescope (JCMT). The observations were 
performed during several periods between 2007 December 27 
and 2008 January 17, with the majority of the data (60% of total 
integration time) being obtained on 2008 January 17. HARP is 
an array of 16 detectors, only 14 of which were operational at the 
time of our observations. The detectors are arranged in a 4 x 4 
grid, separated by 30". The beam-size of each detector is 14" 
(HPBW). The map was obtained in the jiggle-position switch 
mode with the HARPS jiggle pattern, which provides measure- 
ments on a 6"x6" grid. Because of the non-operational detectors 
and to fully cover the interesting area, observations were c allied 
out with different array orientations, mainly with position an- 
gles of PA =45° and PA = 0°. In addition, the map center was 
changed for different observing runs. In consequence, the shape 
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line 


position 


date 


PA" 


spec. res. 




?0N 


rms'' 








["]" 


[km s-'] 


[K] 




[K~ 


'"CO(l-O) 


P32^ 


2008.04.02 


13 


0.081 


406 


lm20s 


1.78 




pio" 


2008.04.15 


58 


0.081 


390 


14m20s 


0.59 


''CO(l-O) 


P32'' 


2008.04.05/06 


13 


0.085 


360 


lh40m 


0.09 




PIO'' 


2008.04.09 


13 


0.085 


390 


IhlOm 


0.13 


HCO+(1-0) 


PIO' 


2008.04.14 


45 


0.105 


240 


54m 


0.09 



" Position angle of BEARS 
* DSB values 

in the 7,^1, (SSB) scale, per the specified resolution element 
a = 07''04"'05^40, 6 = -03°50' 13'.'0 
a = 07''04"'05^52, 6 = -03°50'0a'0 



of the resultant map is irregular and in terms of sensitivity the 
map is rather inhomogeneous. 

The central sky reference (off) position was the same as 
in the NRO observations, i.e. at a = 07''05'"25';70, 5 = 
-04°03'25'.'0. Each jiggle position had a separate off position. 
Observations were performed in good weather conditions with 
a typical atmospheric opacity at 220 GHz of tcso - 0.10 (as 
measured by the tau-meter at the nearby Caltech Submillimeter 
Observatory). The receiver noise temperature (no atmospheric 
contribution) was in the range 109-201 K for different detectors. 
Pointing was performed at 345 GHz regularly every 1 h on strong 
stellar sources and the pointing accuracy should be within 4'.'5 
(3o"). The calibration was controlled by observations of standard 
targets, e.g. CRL618 and IRC+10216, and its stabiHty was better 
than about 10%. 

As a back-end, we used the ACSIS autocorrelator in the con- 
figuration, which provides a 250MHz (217km s"') bandwidth 
and a 31 kHz (0.027km s ') channel resolution. The data were 
reduced with Starlink software, in particular with the Kernel 
Application Package {KAPPA). Baselines of up to third order 
were subtracted from all spectra and the data were regridded 
onto a data-cube with a 6"x6" spatial sampling. The resultant 
data were then converted to a main beam brightness temperature 
with main beam efficiency of //nib = 0.63. 

2.4. IRAM observations in fhe ^^CO(1-0) and (2-1) iines 

For our analysis and discussion, it is impoitant for us to con- 
sider the data reported in Kamiiiski (2008). These observations 
were obtained in 2006 September 27-28 with the IRAM 30-m 
telescope using single-pixel receivers in the frequency switching 
mode. The data consist of very sensitive measurements in the 
'^CO(l-O) and (2-1) lines for 13 positions within the echo re- 
gion, as shown i n Fig.[T] More te chnical details about the spectra 
can be found in lKamifls kil JIooa) . 

2.5. HERA observations obtained in fhe J=2-1 line of '^CO 
and^^CO 

The last spectroscopic observations were performed on 2009 
January 29 with the IRAM 30-m telescope in the '^CO(2-l) 
(230.54 GHz) and '^CO(2-l) (220.40 GHz) lines The HERA 
(HEterodyne Receiver AiTay) (ISchuster et al.ll2004l) was used in 
the stare mode. This array consists of nine detectors arranged in 
a 3x3 grid with a separation of 24" between neighbouring de- 
tectors. In our observations, the central detector was centerd on 
either PIO or P32 with the position angle of the array being set to 
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"2000 



Fig. 1. Selected regions observed in CO lines are shown on top 
of the optical image of the light echo. The red frame shows 
the main region covered by the CO(3-2) map obtained with 
HARP/JCMT in 2008 January. The green circles represent posi- 
tions observed with the single-pixel receivers on the IRAM 30-m 
telescope in 2006 September These positions are labelled with 
numbers introduced in Kamiiiski (2008). The sizes of the circles 
correspond to beam sizes (HPBW) at the observed frequencies, 
i.e. the larger ones represent the CO(l-O) data, the smaller ones 
represent the CO(2-l) data. The blue and magenta circles corre- 
spond to HERA observation obtained with the IRAM 30-m tele- 
scope in the CO(2-l) line in Jan. 2009. The blue and magenta 
colours correspond to observations with the aiTay centerd at po- 
sitions PIO and P32, respectively. The position of V838 Mon is 
indicated with a star. The background image was obtained with 
HST/ACS in the F814W filter on 2006 September 10. 



zero. The observing mode was PSW with an off position located 
600" east from the observed on position. The integration times, 
system temperatures, and the resulting rms values (in units of 
Tmb) are given in Table [3] 

The VESPA autocorrelator was used as a back-end. The total 
bandwidth was 71 MHz (95 km s"') with a spectral resolution 
of 78 kHz (0.1 km s"'). The beam-width at both observed fre- 
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Table 3. Characteristics of the HERA observations from 29 Jan. 
2009. 



line 


position 




foN 


rms* 






[K] 


[min] 


[mK] 


'^CO 7=2-1 


P32 


580 


44 


118 




PIO 


635 


44 


133 


"CO 7=2-1 


P32 


390 


210 


34 




PIO 


400 


210 


35 



° SSB value 

* An average value for all of the 9 beams per resolution element 



quencies is practically the same with HPBW=H". Data were 
calibrated with the standard chopper-wheel method. The point- 
ing was controlled by regular observations of strong continuum 
sources and should be of higher accuracy than 4'.'5 (3cr). Some 
observations were obtained at very low source elevation (~18°). 

The data were reduced using the CLASS package, which 
is part of the GILDAS software collection. Baselines were 
corrected by subtracting low-order polynomials from indi- 
vidual scans. For each position, the scans were averaged 
and transformed to the Tmb scale using efficiencies from 
Fnium & Mauersbergeii (l2007l) . i.e., /7n,b(220 GHz)=0.59 and 
77n,b(230 GHz)=0.57. 

2.6. HST optical observations of tlie liglit eclio 

Multi-epoch observations of the light echo have been performed 
with HST on several occasion s since 2002 and were partially 
published in Bon d et all (l2003h and iBondl (l2007h . For the sake 
of data presentation and the analysis, we used the data obtained 
on 2006 Sept. 10 with the Advanced Camera for Surveys (ACS). 
The CCD frames were downloaded from the Multimission 
Archive at the Space Telescope Science Institute (MAST). They 
were reduced with the standard pipeline and reprocessed with 
the multidrizzle procedure implemented in the STSDAS package 
of pyRAF. 

3. Results of the observations 

3.1. Maps of tile ^^CO(1-0) line 

A careful analysis of the ''CO(l-O) data-cube uncovers three 
kinematically distinct cloud complexes in the mapped region. 
They can be assigned to the velocity ranges 47.0-52.0 km s"', 
52.0-54.5 km s"', and 55.0-58 km s"'. Integrated intensity maps 
for these velocity intervals, together with a map of the rms noise 
distribution of the cube, are shown in Fig. |2] The maps were 
obtained with a 7'.'5 x 7'.'5 grid spacing. 

The clouds that appear in the velocity range 47.0- 
52.0 km s"' are located in a declination strip of (approximately) 
-3°51'13" < 6 < -3°55'30". The strip is seen along the 
whole map and forms an elongated complex, which probably 
extends outside the mapped region. The '^CO(l-O) emission in 
this molecular region is weak, and has multiple local maxima 
centers with peak temperatures reaching ~3 K. This molecular 
complex is very diffuse and consists of numerous patchy clouds. 

In the velocity range 52.0-54.5 km s~\ only one cloud can 
be seen. It is located north with respect to V838 Mon. Although 
certainly extended, the cloud has a compact structure. It resem- 
bles a comet with a core close to V838 Mon, and has a length 
of 93" (as measured for an isophot at 10% of the intensity peak) 



along its longer axis. It is approximately aligned in the north- 
south direction (see Fig.|6]l. The integrated intensity peak is lo- 
cated at a = 07''04'"05!52, 6 = -03°50'0a.'00. This position is 
52" away from V838 Mon. At a majority of positions, the line 
shapes are Gaussian and have no structure, although some lines 
are slightly asymmetiic. The typical full width at half maximum 
(FWHM) is AV = 0.9 km s ' . 

In the under-sampled part of our map, at approximately a = 
Q-jhQ^in^y^ 5 = -03°51'00", an emission region can be found at 
velocities 55.0-58.0 km s"' . This cloud is seen in the very noisy 
part of the map, close to the map's edge, where it is difficult to 
characterize the feature in more detail. 

3.2. The HARP map in the ^^CO(3-2) line 

The region observed with HARP is much less extended than the 
map obtained in CO(l-O) and corresponds to a central portion 
of the large OTF map. The HARP observations in the '^CO(3-2) 
line reveal an extended region of emission in the velocity range 
52-54.5 km s"' and in approximately the same sky region as the 
cloud found in CO(l-O) data at the same velocity range. A map 
of integrated intensity and a corresponding map of rms noise are 
presented in Fig.[3](see also Fig.|6]l. The lines in the emission re- 
gion are very narrow with a typical FWHM of AV = 0.9 km s ' . 

The general appearance of the cloud is similar to that of 
the analogous cloud found in the CO(l-O) map. The feature 
seen in CO(3-2) clearly has a cometary-like structure with a 
core oiiented towards the position of the star The 10% con- 
tour of integrated intensity has a dimension of ~60" along the 
longer axis, and ~40" in the orthogonal direction. The point 
with the strongest emission is located at a = 07''04"'05H, 
5 = -03°50' 13'.'0, that is 39" from the position of V838 Mon. 

3.3. Measurements of single positions observed with BEARS 

3.3.1. '2CO(1-0) 

The observations in the '^CO(l-O) line at PIO and P32 as cen- 
tral positions were very short and intended to be control obser- 
vations. The line was clearly detected only at the two central 
positions. The line parameters are listed in Table The mea- 
sured intensities of these lines are important for understanding 
the uncertainties in the OTF data (see Sect. l3.4.TT l. 

3.3.2. '^00(1-0) 

The BEARS observations in the '^CO(l-O) line only resulted 
in detections for the central detector, i.e. at PIO and P32. These 
detections are at the level of 5cr. The line parameters are listed 
in Table |4] The central velocities of the observed emission lines 
agree very well with the corresponding velocities for ^^CO at the 
same positions. The '^^CO lines seem to be narrower than their 
'^CO counterparts, but this is uncertain for these noisy spectra. 

3.3.3. HCO+ (1-0) 

In the spectrum acquired in the HCO^(l-O) transition, no line 
is detected and only an upper limit to a potential detection can 
be given. For a 0.105 km s"' channel bin, the limit on the peak 
intensity is rmb=0.27 K (3cr). This value is only half of the peak 
intensity of the '■'CO (1-0) line at the same position. 
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Fig. 2. Upper left and right and lower left: Integrated intensity maps in the '^CO(l-O) line obtained with the BEARS instrument at 
the NRO 45-m telescope. The signal was integrated in the velocity range indicated on each map. Contours are drawn for 30, 45, and 
60% of the maximal integrated signal. Lower right: The noise-level distribution in the '^CO(l-O) data cube (rms values are given 
per 0.1 km s ' bin). Dashed lines are drawn for constant declination and right ascension. Note that the frames are not aligned with 
the axis of the equatorial system. The star marks the position of V838 Mon. 



Table 4. Characteristics of the 7=1-0 lines of '^CO and '^CO 
observed at positions PIO and P32 with the central BEARS re- 
ceptor in Apr 2008. 



position 


Vlsr AV Ico rms 
[kms-'] [kms-'] [K] [Kkms-'] [mK] 




"CO(l-O) (OTP) 


PIO 
P32 


53.2 0.9 6.10° 5.59 289.9^ 

53.3 0.9 6.46° 6.33 287.4~ 




'-CO(l-O) (PSW) 


PIO 
P32 


53.4 0.8 8.13° 6.71 590.0° 

53.5 0.9 8.07° 6.77 1782° 




'-'CO (1-0) (PSW) 


PIO 
P32 


53.3 0.7 0.59" 0.45 129.0^ 
53.1 0.5 0.51* 0.27 93.8" 



° per 0.081 km s ' channel 
per 0.085 km S-' channel 



3.4. Measurements for single positions observed witli HERA 

Spectra acquired with the HERA array in 2009 January are 
shown in Figs.|4]and|5] Line measurements are shown in Tables 
|5]and|6] Among the positions observed with HERA centerd on 
PIO, no emission in '^CO(2-l) was detected at the offset (24,- 
24). At positions (-24,-24) and (-24,24), very weak emission 
is seen with an integrated intensity below 1 K km s-' . At the re- 
maining six positions, the lines are strong, with the strongest line 
observed at PIO. The line profiles analysed at the original reso- 
lution are irregular. In at least two cases, i.e. at PIO and (-24,0), 
the profile consists of two components forming a tight blend. For 
two positions with strong '"CO emission, the '^CO line was also 
detected. While at PIO the line is strong, the line at (24,0) is de- 
tected at the level of 3cr. The emission lines of '^CO are clearly 
narrower than the corresponding lines of '^CO. 

In the case of observations centerd at P32, the '^CO(2-l) 
line was detected only at the six northern positions, i.e. no lines 
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Fig. 3. Left: Integrated intensity map of the CO(3-2) emission in the range 52-54.5 km s ' , as observed with HARP/JCMT. Contours 
are at 5, 10, 25, 40, 55, and 70% of the peak emission. Right: The corresponding rms map, where rms is given per 0.1 km s"' bin. 
Dashed lines are drawn for constant declination and right ascension. Note that the frames are not aligned with the axis of the 
equatorial system. 



Table 5. Measurements for lines observed with HERA/IRAM 
in Jan. 2009 with PIO as the central position. The values of T™" 

are for channels of 0.102 km s '('^CO) and 0.106 km s '('^CO). 
Offsets are given with respect to PIO. 



Offset wrt PIO 


Vlsr 


AV 


mb 


^co 


(Aa",A5") 


[kms-'] 


[kms-i] 


[K] 


[Kkms-'] 






'^CO(2 


-1) 




-24 


-24 


53.3 


0.9 


0.52 


0.49 


-24 





53.3 


1.3 


1.63 


2.33 


-24 


24 


53.5 


1.0 


0.47 


0.51 





-24 


53.1 


1.0 


2.68 


2.77 








53.4 


1.0 


6.10^ 


6.46° 





24 


53.4 


0.8 


2.3 r 


2.01 


24 





53.3 


0.9 


1.23 


1.16 


24 


24 


53.4 


0.8 


1.87 


1.68 






'"CO(2 


-1) 










53.3 


0.6 


0.67 


0.42 


24 





53.2 


0.6 


0.10 


0.06 



" The line profile is not Gaussian. The peak of the profile is 
at T|™=5.47 K and the signal summed over the profile gives 
/co=6.15 Kkms-' 



are seen at (24,-24), (0,-24), and (-24,-24). The observed line 
profiles are again irregular, but the line is clearly double only at 
(0,24). The isotopologue line '^^CO(2-l) is only detected at P32 
and (0,24). These lines are almost half as narrow as their '^CO 
counterparts. 

3.4.1 . Comparison between different data sets for tine 
'^CO(1-0) line and calibration issues 

To help us understand the quality of the data discussed in this 
paper, we compared spectra obtained with different telescopes 
and observing modes. The comparison is made for the '^CO(l- 
0) line, which was observed at a few positions over several ob- 
serving runs. The positions PIO and P32 were observed with the 
BEARS aiTay in the OTP and PSW modes. Two positions with 
strong lines at offsets (-8,-14) and (-8,10) with respect to P32 
(named respectively Off3 and OfF7 in Kaminski (2008)) were 
observed in the frequency-switching mode with single pixel re- 
ceivers at the IRAM 30-m telescope, as well as with the NRO 



Table 6. The same as in Table |5] but where P32 is the central 
position. 







'-C0(2 


-1) 




Offset wrt P32 


Vlsr 


AV 


mb 




(Aa",A^") 


[km s-'] 


[kms-'] 


[K] 


[Kkms-'] 


-24 


53.3 


1.2 


0.81 


1.04 


-24 24 


53.3 


0.8 


0.55 


0.48 





53.1 


1.0 


8.38 


8.54 


24 


53.5 


1.0 


4.67 


4.74 


24 


53.3 


0.8 


1.21 


1.02 


24 24 


53.4 


0.8 


2.76 


2.47 






'^C0(2 


-1) 







53.1 


0.6 


0.85 


0.56 


24 


53.5 


0.5 


0.26 


0.13 



45-m telescope in the OTP mode. Spectra for these four positions 
were extracted from the gridded OTP data cube. Differences in 
coordinates of the extracted OTP spectra with respect to those 
defining the four positions should be much smaller than the for- 
mal pointing errors of the OTP observations. 

We first compare the data obtained with BEARS in the PSW 
and OTP modes. The emission lines are on average 1.4 times 
weaker in the OTP data than in the PSW observations. This dif- 
ference exceeds the formal calibration eiTor expected for our 
BEARS observations (~14%). Some velocity shift between the 
two data sets can also be seen. 

When the OTP spectra are compared with the IRAM data, 
the discrepancy is even more pronounced. The OTP line inten- 
sities are 1.6-1.9 weaker than in the IRAM data. Moreover, the 
line profiles for the Off3 position are significantly different in 
the two spectra. Although many factors may be responsible for 
the noticeable differences, the problematic pointing of the 45- 
m telescope may be the main source of errors. This problem 
may be applicable to both OTP and PSW observations and all 
the BEARS results should be treated with special care as it is 
probable that they are less accurate than the formal uncertainties 
indicate. 

3.5. ^^CO(3-2) spectra extracted for selected positions 

In the analysis of line intensities at positions PIO and P32, we 
extracted relevant spectra from the HARP/JCMT data cube. The 
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45 50 55 



60 



Vlsr [km s" ] 

Fig. 4. Spectra acquired using the HERA array in Jan. 2009 with PIO as the central position. Spectra obtained in '^CO(2-l) are 
presented with a blue line, and those for '^CO(2-l) are shown with red line. The offsets of the observed positions are indicated in 
the left upper corner of each panel and are given with respect to PIO. 



i4 




45 50 55 60 
Vlsr [km s"'] 

Fig. 5. The same as in Fig.|4] but where P32 is the central position. The offsets are given with respect to P32. 



differences between the coordinates of the observed positions 4. Analysis: physical parameters Of the molecular 

and those defining PIO and P32 are smaller than the pointing er- clOUd 

rors at JCMT. The results of line measurements for these spectra 

are given in Table [T] 4.1. Distribution 

The rich obser vational material presented in Sect. |2] and 
iKam iiiski' (2008*) enables us to describe quite accurately the spa- 
tial distribution of molecular gas seen in the area of the sky cov- 
ered by the light echo. The maps made in the CO(l-O) and (3-2) 
transitions show a compact region of molecular emission tens of 
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Table 7. Characteristics of the '^CO(3-2) emission lines in spec- 
tra extracted from the data cube obtained with HARP/JCMT. 

are given for a bin of 0.1 km s"^ 



The values of rms and T™" 

mb 



position 


Vlsr 


AV 


mb 


^co 


rms 




[kms-'] 


[kms-'] 


[K] 


[Kkms-'] 


[mK] 


PIO 


53.4 


0.9 


3.58 


3.44 


217 


P32 


53.2 


1.0 


5.29 


5.60 


219 



arcsec north of V838 Mon. The cloud seen in both transitions is 
shown in a contour map of Fig. |6] The observations used to pro- 
duce the maps are rather limited and show only regions of main 
emission in CO. The more sensitive single-pointing observations 
provide an opportunity to constrain even more tightly the gas dis- 
tribution on the sky. All the positions where point observations 
led to a detection of the CO(l-O) or CO(2-l) lines are marked in 
Fig.|6] We first analyse the single point observations in which the 
CO(l-O) line detected, which are indicated in Fig. |6] with purple 
circles. These are the IRAM observations reported in iKaminskH 
(l2008h . Four of those positions are located outside the region 
of main molecular emission, which is defined by the contours 
of the OTF data. The emission found at the position labelled as 
OffS has a central radial velocity of ~48 km s ' and is, most 
probably, not related to the echo material (see Sect. 15. H . The 
emission lines detected at the three remaining positions outside 
the contours (Off2, Off4, star position) are weak (r,nb <0.12 K) 
and reside at a radial velocity of 53.3 km s"', the same as for 
the main emission region in the OTF data. Although these posi- 
tions were covered by the OTF observations, the emission lines 
therein are too weak to be seen in the noisy OTF spectra. The 
deep IRAM integrations prove however that the region of CO(l- 
0) emission is more extended than shown by the OTF map. In 
particular, some weak emission is seen very close to the position 
ofV838 Mon. 

The positions with detections of the CO(2-l) line are shown 
in Fig. |6]with green pluses. The emission at Off8 is not of in- 
terest for the same reasons as given above. All the positions 
with CO(2-l) detected, within the beam-sizes (c.f. Fig. [T]), re- 
side within the contours of the CO(l-O) emission. Even though 
long integrations were obtained at several positions outside the 
main cloud region, no 7=2-1 emission was found outside the 
contours. In particular, there is no emission at the three south- 
em positions where weak lines of CO(l-O) were found (01i2, 
Ofr4, position of V838 Mon). Remarkably, the IRAM observa- 
tions obtained at the position of V838 Mon with an integration 
time of 6 h gave only an upper limit to the CO(2-l) emission 
of 3cr=13.9 mK km s This means that the CO(2-l) line is at 
least 10 times weaker than the corresponding CO(l-O) line at 
the same position, and 20 times weaker than the CO(l-O) line 
at the other two neighbouring positions (OfF2 and Off4). It may 
be concluded that the region in the vicinity of V838 Mon lacks 
CO(2-l) emission. 

4.2. Modelling the physical parameters of the cloud 

Our observational material provides measurements of multiple 
transitions at several positions located within the echo region. 
At least three CO lines were measured for the four positions 
PIO (5), P32 (5), Off3 (3), and Off? (3) (where the numbers in 
brackets indicate the number of measured lines). These measure- 
ments can be used to model the physical conditions within the 
CO cloud. We compare the measured intensities with those sim- 
ulated for a model cloud with physical conditions varying over 
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Fig. 6. Contours of '^CO emission overplotted on the light echo 
image from HST (obtained in Sept. 2006, filter F814W). The 
blue contours show the '^CO(l-O) emission integrated in the 
velocity range of 52-54.5km s"' (0.99, 1.35, 2.01, 3.26, and 
5.66 K km s '), while cyan contours correspond to '^CO emis- 
sion integrated in the range 47-51 km s"' (1.19, 1.34, 1.57, 1.75, 
2.07, 2.40, and 2.87 K km s"'). The red contours show emis- 
sion in the '^CO(3-2) line integrated over 52-54.5 km s"' (0.31, 
0.51, 1.07, 2.10, and 3.88 K km s"'). The position of V838 Mon 
is indicated by the star. The green '-h' symbols mark positions 
with detected emission in the '^CO(2-l) line in the I RAM sin- 
gle position observations (Off'3, Off7, and OfF8, as in iKamiriskil 
(2008)). The 'x' symbols mark positions with the strongest 
emission from the '^00(1-0) line (blue, PIO) and ^^CO{3-2) 
line (red, P32) integrated in the 52-54.5 km s"' range. 



a broad range. For rad iative transfer calculatio ns and the simula- 
tions, we used Radex (van der Tak et al. '2007) (in a version from 
2008 Nov. 20). The program and the input data are described in 
Sect. 14.2.11 The best-fit set of simulated line intensities was de- 
cided using the test, where we used the formula 



X 



/ robs 
^CO,/ 



rmod \ 

^CO. 



C0,( 



(1) 



The summation is made for expressions corresponding to each 
transition measured, where I™^. is the simulated intensity, /^q . 

is the measured intensity, and A/^q . is the measurement uncer- 
tainty. As the uncertainty, we took the geometric mean of the 
absolute calibration error and the value of 3cr().iAych, where cro.i 
is the rms of the spectrum smoothed to the channel resolution of 
AVch=0.1 km s"'. The absolute calibration errors were assumed 
to be 25% for BEARS/NRO, 20% for HARP/JCMT, 15% for all 
the IRAM data at 230 GHz, and 10% for all the IRAM CO(l-O) 
data. The assumed BEARS errors are higher than the previously 
indicated formal calibration errors. This is related to the calibra- 
tion inconsistency discussed Sect. l3.4.TI 

Although in our calculations we included all the positions 
with multiple CO lines measured, we discuss here only results 
for positions PIO and P32, where five transitions were observed. 
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and therefore the constraint on physical parameters are strongest. 
To make the analysis more independent of the uncertain calibra- 
tion of the BEARS measurements (see Sect. 13. 4. Tb . we compared 
observations and simulations (/) for all the measured lines, (//) 
by excluding the BEARS data for '^CO(l-O), and (Hi) indepen- 
dently of all the BEARS data. 

4.2.1. The code 

Radex is a non-LTE radiative transfer code, which uses the es- 
cape probability formulation. We used the program in the uni- 
form sphere mode, in which the medium is treated as static, 
spherically symmetric and homogenous. 

Input spectroscopic and collisional da ta required by Radex 
were taken from the LAMDA database (ISchoier e t al. 2005^. 
where collisional rates for CO-H2 are taken from .Yang et al] 
( I2OIOI) . Though Radex enables computations to be made for 
many collisional partners, we considered only collisions with 
H2. However, we note that in our case atomic hydrogen can also 
be an important collisional partner of CO, but that almost noth- 
ing is known about the H/H2 ratio in the cloud. We also ignore 
collisions with He, for which LAMDA does not provide collision 
rates. The LAMDA database provides CO-H2 collision rates for 
ortho- and para-H2 molecules. In our computations, the ortho- to 
para-H2 ratio was assumed to be thermalized at the given kinetic 
temperature of the model. As the background radiation field, we 
implemented only the cosmic microwave background with the 
black-body temperature of T\,„ - 2.725 K. 

The entire radiative transfer in Radex is performed on rect- 
angular line shapes, i.e. the optical depth is not changed over 
the profile. To compare the computed line intensities with the 
measured velocity-integrated line intensities, the computed radi- 
ation temperature, Tr, is corrected within the program by a fac- 
tor ( ^Jn|2 Vin2) x AV = 1 .064 AV, where AV is the user-defined 
Gaussian half-width of the line. The typical observed widths for 
the '^CO lines are AV - 0.9 km s ', while for the '^CO Unes 
AV - 0.6 km s ' (cf. Sect.O. We interpret the broadening of the 
^^CO lines as the effect of their moderate or high optical thick- 
ness (t > 1). The radiative transfer calculations within Radex 
were performed with AV - 0.6 km s ' for all lines, but the final 
integrated intensities were calculated with tsV - 0.9 km s"' for 
the ''CO lines and AV = 0.6 km s ' for the '^CO fines. 



4.2.2. The models and results 

Using Radex, we generated several grids of models with wide 
ranges of physical parameters. The initial grid was constructed to 
place strong constraints on the CO column densities, A^('^CO). 
It was calculated for column densities in the range lO'"* - 



10 cm with values changed in a step of 0.25 dex (deci- 
mal exponent). For each value of CO column density, simu- 
lations were performed for H2 densities from the set of val- 
ues given by n(H2)=[2, 4, 6, 8, 10]xlO* cm"^ where k =2, 
3, 4,..., 8, and for kinetic temperatures from the range 7'idn=5- 
200 K separated by a step of 5 K. For the '"'CO lines, the col- 
umn density was assu med to be given by th e standard ratio 
A^('^^CO)=A^('2CO)/60 jLanger & Penziasl[l99"3 ). For this grid, 
the minimum of is consistently reached for both positions at 
log[A?('2cO)/cm-2]=16.25 or log[A^('2cO)/cm-2]=16.75 when 
all the BEARS data are excluded. 

As discussed later in Sect. 14.31 the analysed molecular gas 
is most probably located in the outer parts of the Galaxy, 
where the isotopologue ratio A^('^CO)/A^('"'CO) is known to 




3.5 4 4.5 5 5.5 

log[n(H2)/(cm^)] 




log(Z^) 



1.5 



0.5 



3.5 4 4.5 5 5.5 6 6.5 

log[n(H2)/(em-3)] 

Fig. 7. Distributions of log(x^) in the [Tkm, n(H2)]-plane for the 
4 measurements at PIO (top) and P32 (bottom) and '^CO column 
density of log A^= 16. 25. The blue contours are drawn for x~-^-^ 
and 2.5, while the red contours mark;t'^=LO and L2. The pluses 
indicate the minima. 



be higher than the standard value in the solar neighbour- 
hood (e.g. iBrand & Wouterloo^ Il995l) . We therefore gener- 
ated a grid of models that is similar to the one described 
above but with NC^CO)^NC^CO)inO, where the high iso- 
topologue ratio was taken after Brand & Wouterloot (1995)Q 
This grid suggests for both positions that the column density is 
log[A^('^CO)/cm"^] = 16.5, but the fit to observations is slightly 
worse than in the previous grid of models. 

All subsequent simulations were performed with 
A^(13CO)=A^('2CO)/60 and log[A?(i2cO)/cm-2] = 16.25, as 
these values result in best fits to the observations. All the initial 
models gavsx^ minima at rkin=15±5 K. 

With A^('^CO) and the CO isotopomer ratio fixed, we gen- 
erated a denser grid of models at 2 K< <80 K in steps 
of 0.5 K, and densities given by n(H2)=r ■ 10* cm"-', where 
r=l,L5,2,. . .,9.5 and k-2,3,. . .,1 . The simulated line intensi- 
ties were compared with the measured intensities excluding 
the BEARS data for '^00(1-0). The distribution of in the 
density-temperature plane is shown in Fig. |7] 

The minimum of the x^ distribution for PIO is at x^-^-^ 
and indicates that «(H2)=9.5-10^ cm~^ and rkin=12 K. For 
P32, we get a broader but deeper minimum at x^-^-'^^ giving 



' The high A'('^CO)/A'('''CO) ratio does not necessarily have to be 
related to the location of the cloud in the outer Galaxy, but may be 
related to the isotope-selective nature of the phot o-dissociation proc ess 
of the cloud by the interstellar radiation field (see lVisser et alj2009h . 
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n(H2)=3-10^ cm"-' and rkin=14 K. Taking into account uncer- 
tainties mainly in the derivations of the column density, the phys- 
ical conditions towards PIO and P32 may be considered to be 
identical and represented by the values of n(H2) ~ lO'* cm"^ and 
rkin~13 K. In these conditions, the '^CO lines arising from levels 
7<3 are of moderate thickness, i.e. 1< r <10, and some broad- 
ening with respect to the optically thin '^CO lines is expected. 
This justifies our modelling procedure described in Sect. 14.2.11 

For the derived physical parameters, we also calculated the 
intensities of the HCO^(l-O) line assuming the molecular abun- 
dance of HCO+/H2=2-10-'' (Turner 2000). These predicted line 
intensities are of the order of 10"^ K, so far below our detec- 
tion limits. The observed cloud is not dense enough to produce 
HCO^(l-O) line fluxes, which would be easily observed. 

4.3. Kinematical distance and location in tiie Galaxy 

The molecular cloud of interest has a very well defined radial 
velocity. By fitting a Gaussian to the profile of the emission 
averaged over the whole cloud, we get a central velocity of 
Vlsr=53.34 km s"' . The cloud has Galactic coordinates /=217°8 
and b=l°.Q, located close to the outer plane of the Galactic disk. 
Using rotational curves for this part of the Galaxy, we can es- 
timate the ki nematical dista nce of the cloud. For the rotation 
curves from [Brand & Bhtzl fl993l) and iFich et all (119891) . we 
found that d-b.5-1 kpc. Possible streaming motions make this 
result uncert ain by about 1 kpc. F or the fully empirical veloc- 
ity field from lBrand & Bhtzl (1 19931) . we get c/=7.0±0.3 kpc. The 
derived values suggest the cloud is located in the Outer Arm 
(Norma-Cygnus Am; d ^6.5 kpc) (fRusseill 120031) . The heHo- 
centric distance of 7 kpc corresponds to the Galactocentric dis- 
tance of 14 kpc, which indicates that the cloud is located on the 
outskirts of the Galactic molecular disk. 

4.4. Column density of H2 and Ay 

To classify the molecular cloud and reliably constrain its relation 
to the light echo material, an estimate of the total column density 
of H2, or equivalently the reddening Ay, is required. The results 
of Radex modelling can be used for this purpose. For the consid- 
ered lines of sight, we obtained log[A^(CO)/cm"^] = 16. 25+0.25. 
Assumin g the standard CO a bundance (with respect to H2) of 
~ 10"* dWilson et al.1 12OO9I) . this result can be converted to 
log[A^(H2)/cm"^]=20.25+0.25. The error does not take into ac- 
count the uncertainty in the CO abundance, which can vary by a 
factor of about 2. 

The total column density of mole cules can be derived inde- 
pendently using the X-factor method (Blo emen et al.|[T986h . For 
an integrated intensity of the CO(l-O) line, Ico, the hydrogen 
column can be estimated from 

A^(H2) = Xco ■ J r,„b(''CO) dV = Xco/co. (2) 

We take ^co=2.8 (in units of 10^° cm^^ K"' km"' s; 
iBloemen et alJ l 1986), although values as high as 6 can be found 
in the literature for the outer Galaxy (see e.g. Kaminski et al. 
I2007h . For positions PIO and P32, we get A^(H2)=(1.6 - 1.8) ■ 
10^' cm"^, which should be representative of the densest re- 
gions of the cloud. This result is one order of magnitude higher 
than for the results of Radex modelling, but the actual errors 
in both estimations are comparable to this difference. For the 
emission profile averaged over the entire cloud (i.e. within the 
isophot at 10% of the peak in the OTF data), we measured 



(^co)=2.88 K km s" , which gives an ,, average" H2 column den- 
sity of 8.1 ■ 10^" cm"^. We assume that this value most closely 
represent the whole molecular cloud. Using the s tandard con- 
versio n formula, N{\{)IEh-v -5 .8 10^' cm"^ mag ' (iBohlin et al.l 
[T9781 iRachford et al.l I200I . we can calculate the correspond- 
ing reddening Eb~v- Ignoring the atomic contribution to the to- 
tal hydrogen column, i.e. A?(H)=A^(HI)+2A^(H2) ~ 2A^(H2), we 
obtain £'b-v=0.28 mag. For the standard ratio of visual to total 
extinction of Ry =3.1, this corresponds to a visual extinction of 
Ay-Q.%1 mag. For the column densities at positions PIO and P32 
estimated with the X-factor method, we get Ay^l.O mag. This 
result shows that even in the densest regions of the cloud the ex- 
tinction is low. This is consistent with the optical HST images 
of the echo region, which show distant galaxies shining through 
the cloud (e.g. at or = 07''04"'05^;2, 6 = -03°49'40"). 

4.5. Mass of the cloud 

The mass of the cloud can be estimated using our '^CO(l-O) 
data. First, we estimate its virial mass. For a cloud with a spatial 
density profile of the type p(r) cc r^^, with radius R (expressed 
in pc), and velocity dispersion AV (in km s~'), the virial mass 
(expressed in Mq) is given as dMacLaren et al.|[l988l) 

Mvir = 126/?Ayl (3) 

For the 10% isophot of the CO(l-O) emission and the distance 
of d-6.5 kpc, the radius is R-\.5 pc. The emission profile av- 
eraged over the entire cloud has a width of Ay=0.9 km s 
These values give a virial mass ofMvir=153 Mq. For a flat den- 
sity profile, i.e. p(r)=const, which is less likely, the result would 
be Mvir=255 M©. 

An alternative way of deriving the mass of the cloud is by 
using the X-factor method. The mass of molecular matter within 
the cloud can be expressed as 

M(H2) = Xco^d^mH,{Ico), (4) 

where Q is the solid angle of the emission region, d is the dis- 
tance to the cloud, toh, is the mass of the H2 molecule, and (/co) 
is the average CO(l-O) line intensity within the cloud. For the 
solid angle expressed in arcmin^, distance in kpc, line intensity 
in K km s"', and Xco=2.8-102" cm^^ K"' km"' s (.Bohlin et all 
1 19781) . Eq. (|4]i can be rewritten as 

M(H2) = 0.4Qflf2(/co>Mo. (5) 

For the solid angle Q.-nr^ -1.9 arcmin^, the distance of 6.5 kpc, 
and average intensity of (/co)=2.88 K km s"', we obtain a mass 
of M(H2)=92 Mq. Correcting this value for the presence of 
He, which constitutes about 20% of the hydrogen mass, we get 
M(H2 + He)=l 1 1 M0. Owing to a gradient in the value of Xqo 
towards the outer Galaxy, the mass may even be twice as large. 
Taking into account all uncertainties, the derived mass is consis- 
tent with the virial masses found above. 

4.6. Classifications of the cloud 

The molecular cloud seen within the light echo region has all 
the characteristics of an interstellar molecular cloud. In an at- 
tempt to classify it into one of the kno wn groups of interstellar 
clouds, we refer to the classification in ISnow & McCaiil (l2006l) 
and van Dishoeck et al. ( 1993). 

The densities and column densities (or Ay) of molecular mat- 
ter in the observed cloud are too low to classify it as a dark cloud. 
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The maximal density of our cloud of ~ 10^ cm"^ is in the range of 
average densities typical of dark clouds. On the other hand, the 
derived density of the cloud is too high to consider it as a diffuse 
cloud, for which densities in the range n(H2)=100-500 cm"^ are 
expected. In addition, the derived temperature of ~13 K is too 
low compared to the range 30-100 K that is typical of diffuse 
clouds. We therefore classify the cloud as translucent, which 
means in-between being dark and diffuse. The derived extinc- 
tion of Ay ~1.0 mag, which does not take into account a contri- 
bution of atomic matter, is consistent with the range ~l-2 mag 
characterising translucent clouds. The derived mass, dimension, 
densities, and temperatures also support this classification. 

Observations of molecular clouds show that they have a n 
onion-Uke structure (ISnow & McCallll2006l:IWolfire et al.ll2010h . 
Translucent clouds are surrounded by a layer of gas with physi- 
cal parameters typical of diffuse clouds. A layer beyond this may 
exist, which is dominated by very diffuse atomic gas. It seems 
that the molecular cloud observed in the field of the light echo 
has such a complex structure. The CO emission detected close 
to (positions Off2 and Off4) and at the stellar position, i.e. out- 
side the main region of molecular emission, corresponds to these 
outer diffuse regions. A calculations analogous to that described 
in Sect. l4.4l vields Ay ~0.03 mag for this region (the contribution 
of an atomic component to Ay has not been included). 

5. The CO cloud and the light echo material 

The observed molecular emission centerd at 53.3 km s"' coin- 
cides spatially with the light echo seen in the optical images. 
This does not necessarily mean that the two environments are 
physically related, i.e. they do not have to reside at the same dis- 
tances. Below we argue that the molecular gas and the scattering 
dust are actually manifestations of the same interstellar cloud 
located in the vicinity of V838 Mon and its open cluster 



material. A match between radial velocities of objects along the 
same line of sight is usually taken as strong evidence of the phys- 
ical association of the objects. We therefore postulate that the 
molecular cloud seen at Vlsr - 53.3 km s"' and the dust seen in 
optical echo images belong to the same interstellar cloud. 

We note that the molecular emission found at 
Vlsr ~48 km s"' a the radial velocity close to the maser 
velocity (echo material velocity), but that its spatial distribution 
(see e.g. Figs. |2] and |6]l does not indicate any relation with the 
echo material. This cloud is probably a foreground cloud with 
respect to the scattering medium around V838 Mon. 

5.2. Distance 

The distance to the scattering cloud of dust seen in the echo 
images is very well constrained. The most accurate value was 
found with the polarimetric meth od of ISparksl (1 1 9941) . i.e. 
ii=6. 1+0.6 kpc (S parks et al.l l2008l) . An independent estimate 
comes from fitting the main s equence to the photo metry of the 
V838 Mon's cluster members (lAf^ar & Bondl2007h . which gives 
ii=6.2±1.2 kpc. These values agree within uncertainties with the 
kinematical distance derived for the CO cloud of d-6.5-1 kpc 
(see Sect. 14.3b . This provides further evidence that the dust and 
gas we observe are components of the same interstellai^ cloud. 

5.3. The spatial distribution of molecular gas and dust 

If the dust seen in the light echo and the molecular radio emis- 
sion originate in the same interstellar cloud, why do they dis- 
play such different spatial distributions in Fig. |6]? This question 
is related to the issues of relative distributions of gas and dust 
in molecular clouds and the complex geometry of the light echo 
phenomenon. Both issues are discussed below. 



5. 1 . Radial velocity 

The radial velocity of the light-echo material has been deter- 
mined directly. To our knowledge, there are no spectral observa- 
tions with a quality and a resolution good enough to obtain these 
measurements. There are however some observational findings 
that enable us to put constrains on the radial velocity. The dust 
seen in the optical light echo images is certainly locate d close 
to V8 38 Mon, i.e. within several pc (Tylenda 2004; Sparks et alj 
l2008h . According to our current unders tanding of V838 Mon, it 
is a young object dAf^ar & Bondll2007 l) and the star is expected 
to follow the velocity field of its local interstellar environment. 
If this were the case, the light echo material should have a veloc- 
ity close to the stellar systemic velocity, i.e. yLSR=54+l km s"'. 
This velocity corresponds to the measured central velo city of the 
SiO maser in V838 Mon ( Deguchi et al.ll2005L 12009 '). and was 
confirmed as t he vel ocity of the photosphere of V838 Mon by 
iTvlenda et alJ (120111) . We note that the star is known to be in a 
binary system with a B3V star, but the se paration between th e 
components must be large (A ~250 AU, iTvlenda et al.l l2009h . 
hence the orbital velocity is expected to be small (a few km s"'). 
We conclude that the scattering material of the echo has a radial 
velocity close to the radial velocity of the SiO maser, with an 
uncertainty of about a few km s"' related to the orbital motion 
ofV838 Mon. 

The radial velocity of the observed CO emission appearing 
within the echo region is, within the uncertainties, identical to 
the maser velocity. With respect that written above, this means 
that the molecular cloud has the same radial velocity as the echo 



5.3.1. On the correlation between the spatial distribution of 
dust and molecular gas 

It is widely known that gas column densities are spatially corre- 
lated very well with extinction in molecular clouds. For instance, 
for dark molecular clouds a very good correlation was found 
for 2<Av<30 mag, dow n to spatial scales of 0.1 pc (lAlves et al.l 
ll999HLada et al.lll99"9l) . It is often assumed that the correlation 
continues to the more diffuse clouds with Ay <2. We note that 
the correlation of the two distri butions is also ex pressed by the 
classical relation A^(H)oc Eg^y (iBohlin et al.lll97 8). Remarkably, 
the relation refers to the total hydrogen column density, i.e. 
molecular and atomic gas. A molecular cloud is expected to pos- 
sess an extended 'envelope' of atomic gas of low molecular con- 
tent and a non-negligible content of dust. 

An illustrative example of the relative distribution of dust, 
atomic matter, and molecular gas comes from an analysis of the 
interstellar medium within the Pleiades. The cluster is known 
to reside within a dusty medium clearly seen at ultraviolet and 
optical wavelengths. I n the direction of the young cluster, CO 
emission is observed (i Federman & WillsonI 119841) . as well as 
atomic gas traced mainly by HI emission at 21 cm (Gibso^ 
2007). Comprehensive observations of the cluster and field stars 
show that t he observed d ust and gas are indeed located within 
the cluster (lBregerill987l) . The medium has a low density and 
the derived extinction is in the range 1.0< Ay < 1.6 mag, which 
makes it similar to the translucent cloud identified in the echo 
region. The relative distribution of atomic gas, molecular mat- 
ter, and dust within the Pleiades is compared on the web page of 



12 



Kaminski: CO emission in tiie ligiit echo 



S. Gibsoifl The molecular emission of CO is limited to a very 
small region and appears close to the highest concentration of 
dust. The dust is seen in a much larger area of the sky than the 
CO emission and its distribution correlates very well with the 
distribution of atomic gas. Although the interstellar medium in 
the Pleiades is somewhat extraordinary (the cloud and the cluster 
do not have a common origin and are thought to be colliding), 
it may be a good point of reference for what one can expect for 
the cloud identified close to V838 Mon. The dust can be ob- 
served outside the main molecular region, although the highest 
dust concentration is expected in the main region of molecular 
emission. 

The optical images of the echo show a bright nebula cov- 
ering a large region of the sky, which at late epochs extended 
to 60"-90" (2-3 pc) from V838 Mon. In particular, the echo is 
seen in the region of the main molecular emission (see Fig. |6|l, 
but certainly extends beyond it. In light of what was written 
above, this does not rule out a physical association of the two 
media and is expected for a dusty molecular cloud with an ex- 
tended atomic envelope. The echo images obtained after 2003 
(Bond 2007) show a bright echo region close to our position P32, 
which can be identified with the direction towards the center of 
the molecular cloud. This bright clump seems to be the center 
of an extended optical feature, which has the form of a swirl in 
the late-epoch echo images. The collocation of the center of the 
feature and the center of our molecular cloud is probably not a 
coincidence and reflects the physical associations of the scatter- 
ing dust with the molecular gas. 

5.3.2. A comparison between the densities of tine eclio 
material and tine molecular cloud content 

When the distance to a source is known, calibrated images of a 
light echo in general enable one to reconstruct the dust distribu- 
tion around the source. For a sequence of images, even a three 
dimensional distribution of the dusty cloud can be reproduced. 
Such a project is pending for the light echo of V838 Mon and its 
result should help us to understand the connection between the 
dust and the molecular matter In the present discussion, we esti- 
mate the density of matter for only one chosen region of the echo 
and limit the analysis to one image from 2006 September 10 ob- 
tained with ACS/HST in the F814W filter The local density for 
a Ught echo element can be estimated from 



(6) 



where r is the distance of the scattering element from V838 Mon, 
Bsca is the surface brightness of the light echo element, d is the 
distance to the star, Az is the geometrical thickness of the scatter- 
ing element along the line of sight (related to the duration of the 
illuminating flash), (F*) is the flux density of the illuminating 
source averaged over the entire flash, Q^^a is the scattering effi- 
ciency; (Tg is the geometric cross-section of a dust grain, O is a 
(normalized) phase function for the considered scattering angle, 
0, characterizing the echo element, and /=«gi /«H is the number 
density ratio of dust and hydrogen. The integral in Eq. (|6) rep- 
resents the whole distribution of grain sizes a. For our analysis, 
we chose a bright clump in the echo coincident with position 
P32. This region is located at an angular distance of 37'.'5 from 
V838 Mon, or equivalently at a projected distance of p=l.l pc 
(at the distance d=6 kpc). Surface brightness was measured with 



a circular aperture of radius 1'.'8. The aperture is small enough 
to avoid contamination by field stars. The surface brightness 
was measured on the calibrated image (see Sect. 12.6) in GAIA, 
giving Bsca(F814W) = 1.5 • 10~2' erg cm^^ ^-i ^-i arcsec"^ 
The conversion from counts (electrons) to flux units was done 
with the scaling f actor photlam-1 .Q3-\Q'^^^ ergcm"^ A(e")"' 
(ISirianni e t al. 2005). This measured brightness should be com- 
pared to the source intensity in a band of the same throughput 
as for F814W (/teff(F814W)=8332 A). Most of the photomet- 
ric measurements of V838 Mon during the outburst were ob- 
tained in the Ic band (/ieff(/c)=7869 A). The average brightness 
of V838 Mon during the main flash was Iq ~5.5 mag, which 
can be convertecOto m(F814W)=/c + 1.21=6.7 mag. Using the 
flux-density scale c alibrated on Vega (zero point at 25.5 mag, 
ISirianni et ani2005l) . fliis gives a flux density, <F*(F814W)), of 
1.3 ■ 10"" ergcm-2 s"' A"'. 

The rapid rise in V838 Mon's brightness in Ic started at the 
beginning of February 2002 and the o bject remained very bright 
in this band for about 80 days (see e.g. lGoranskii et al.l2007 ). By 
solving the echo equation, z - I2ct — ctjl, for the moments of 
the beginning and the end of the flash, we get Az=0.05 pc. For 
the moment corresponding to the middle of the flash and using 
the relation r - (p^ + z^Y^^, we get the distance of the scattering 
element from the star of r=l.l pc. From another echo relation, 
9 - arcsin(p/r), we find the scattering angle 0=1O3°4. The local 
density of matter can be now expressed as 



= l.M0-22[Jg,,3(^,ff,fl)o-gO(6i,^ff,fl)/(fl)dfl]"'cm-isr-i (7) 
= l.M0-22 47r[5(ieff)] 'cm-3 = \3-\Q-^\S{A^«)\^cm-\ (8) 



Values of the integral S {A, 6) were calculated in ISugermanI 
("2003^, where the author does not use the normalized form of 
the phase function, hence we make the transformation from 
Eq. O to Eq. ([8]l. We limit our considerations to silicate 
dust with a standard ratio of the selective-to-total extinction of 
Ry — 3.1 and grain sizes, a, in the range 5-1 - 1 yum (cf. 
IWeingartner & Drainell200ll) . ISugermanI (l2003l) presents values 
for the integral 5 for three scattering angles, 0=0°, 30°, and 180°, 
and we assume his results for 0=180° to be representative of 
our case (the phase function changes very slowly between the 
derived, i.e. 0=1O3°4, and the adopted angle). For the value of 
S =210"^^, we then obtain the local density of the echo material 
towards P32, ne - 65 cm"^. 

The density of the echo material calculated form the op- 
tical image can now be compared to the density derived for 
the same line of sight from radio observations of CO lines. 
By modelling the line intensities at P32, we found the density 
of «H ~ 2«(H2)=6-10'* cm-3 (Sect. l4T2l i. The density of the 
molecular cloud is almost three orders of magnitude higher than 
the value derived from the optical image. Thus, it appears that 
the dust seen in the echo image from 2006 is not in the dense 
molecular region of the cloud but in its outer (mostly atomic) 
layers. It seems that at least up until 2006, the inner, molecu- 
lar parts of the cloud had not been penetrated by the light-echo 
paraboloids. Taking into account the geometry of a light echo re- 
vealing mostly foreground material, the main molecular region 
should be located „behind" V838 Mon. 



^ http://www.naic.edu/~gibson/pleiades/vla/comparisons.html 



^ http://www.stsci.edu/hst/wfpc2/analysis/wfpc2_cookbook.html 
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5.3.3. The glow 

An analysis of the HST images of the echo obtained after 2007 
in the F606W and F814W filters reveals the presence of a re- 
gion close to the position of V838 Mon, which is as bright as 
the clump that coincides with P32. This bright patch spreads out 
to ~15" north and 40" south from the star Hereafter, we call 
this region the glow. The flux density of the brightest clumps in 
this region is 1.3 times higher in F814W than for the clump at 
P32. At shorter wavelengths, the region at P32 is the brighter 
one. The red color of the glow is essential to understanding its 
origin. After the main eruption, which ended with a ~4.5 mag 
drop in Ic brightness, V838 Mon has remained fairly bright in 
the Ic band during its entire subsequent evolution. While in the 
V band the object faded to about 15 mag, it was still brighter 
than ~10.5 mag in Ic- About half a year after the eruption, a rise 
in Ic brightness by ~ 1 mag was observed. We propose that the 
red glow observed near V838 Mon does not reflect the light of 
the main eruption, but is instead related to the significant bright- 
ness of the object at red and infrared wavelengths after the main 
outburst. Scattering of light occurs not in a paraboloidal layer, 
but in a ^//eii paraboloid and the corresponding thickness of the 
scattering region, Az, can be large. This situation is analogous to 
an echo observed before the end of a flash, as considered e.g. in 
lTvlendal(i2004.) . 

In the proposed interpretation of the glow, the scattering oc- 
curs at small angles and is therefore very effective, so the emerg- 
ing nebula can compete in brightness with the echo of the main 
outburst observed years after the eruption. It is not possible to 
estimate the local density of the glow without knowing the dust 
distribution within the appropriate paraboloid, but it is likely that 
the effective thickness corresponding to the glow is much larger 
than for the region close to P32. Because of the large thickness 
and high scattering efficiency, the local densities corresponding 
to the glow may be much lower than those derived for P32 in 
the analysis of the optical image. This leads to a conclusion that 
the region close to P32 coincides with the highest density echo 
material seen in the optical. This, in turn, shows that the direc- 
tion toward the center of the molecular cloud coincides with the 
highest density region of the echo, as expected if the two envi- 
ronments are manifestations of the same interstellar cloud. 



6. Discussion: the origin of the echo material 

Observations of the ,, infrared echo" (which we prefer to call the 
thermal echo) with the Spitzer Space Tele scope place initial con- 
straints on the mass of the echoing dust. iBaneriee et al.l (l2006l) 
estimated the mass of the thermally radiating dust to be 1 .6 M0 
for the distance of 8 kpc, which when converted to the distance 
of 6 kpc becomes 0.9 Mq. With the standard gas-to-dust mass 
ratio of 100, the total mass associated with the echoing cloud is 
about 90 M0. The geometry of the thermal echo is analogous to 
the optical echo, so this mass corresponds to a thin paraboloidal 
layer within the cloud. The total mass of the cloud may be higher 
This estimate is however in excellent agreement with our con- 
straints on the cloud mass of 90-150 M0 (Sect. 14.5b . Such a 
massive cloud cannot be the result of mass loss from any sin- 
gle star (or even multiple stars) and is definitely interstellar in 
nature. Below we investigate its origin in more depth. 

V838 Mon is a member of a group of hot stars, which in- 
cludes three stars of spectral types B3-6 and the B3 companion 
of V838 Mon ( Afsar & Bond 2007). AH of these stars are seen 
within the echoing region of the sky but well outside the region 
of the main molecular emission. The angular distances of those 



stars from V838 Mon correspond to distances <1.2 pc, while the 
echo until 2008 was seen at distances reaching at least 3 pc. If 
the cluster is not extraordinarily elongated in the radial direc- 
tion, the stars should reside within the interstellar cloud (but not 
within the molecular part of it). We propose that the dust and gas 
seen in optical and radio observations is material remaining after 
the formation of the cluster 

The cumulative ultraviolet radiation of the group members 
must have an influence on the interstellar cloud. The cometary 
shape of the molecular region seems to support this notion. The 
head of the cometary feature is directed toward the center of the 
cluster, as often observed for clouds exposed t o a stron g radiation 
of ho t stars (self-shadowing eff ect; see e.g. iBertoldi &McKei 
Il990t lLefloch&Lazareffl |1994| ) The cloud is flierefor e in die 
process of a „chemical decay" (iDiaz-Miller et al.ll 19981) . which 
is typical of photon-dominated regions (PDRs). The low column 
density of our cloud corresponding to Ay «1 mag, means that 
it is almost transparent to the interstellar radiation field and the 
radiation of the cluster This implies that it should completely 
decay on a timescale of Myr That we observe the cluster to 
be embedded within the cloud provides strong evidence of its 
youth. Clusters containing stars of intermediate B sub-types are 
indeed not observed in association with m olecular gas if the clus- 
ter is older than 3 Myr (lAllen et al.l2007[ and references therein). 
Moreover, it is predicted that clusters with no massive stars, as 
i n our case, disperse a tomic and molecular gas within 10 Myr 
(iPaUa & Stahlej|2000h . The cluster of V838 Mon and the asso- 
ciated gas should also be in a process of decay in a dynamical 
sense. The relaxation time of open clusters with small numbers 
of stars (i.e. less than 36) is estimated to be short, i.e. a few Myr 
(Allen et al. 2007, and references therein). Since the observed 
cluster is still quite compact (~1.2 pc), its age should not exceed 
a few Myr 

It may be concluded that the observed molecular cloud is 
decaying and the a ge of the cluster embedded within it is a few 
Myr, say 3-10 Mvrl Af^ar & Bond (2007) constrained the age of 
the cluster to be younger than 25 Myr. Our result is consistent 
with their constraint, but suggests that an even younger age is 
more probable. 

The above finding is relevant for understanding V838 Mon. 
Some authors (e.g. .Bond_200 7) have claimed that the matter seen 
in the light echo is circumstellar in nature, which would imply 
that V838 Mon is an evolved object with a long mass-loss his- 
tory. This interpretation contradicts our results. The mass of the 
echo material is too high to be associated with any past mass 
loss events in V838 Mon. Moreover, as outlined above, the as- 
sociation of V838 Mon with the interstellar cloud suggests that 
the star must be a young object. 
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